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ABSTRACT 

We report the detection and monitoring of transient substructures in the radiation-driven winds of five mas- 
sive, hot stars in different evolutionary stages. Clumping in the winds of these stars shows up as variable, narrow 
subpeaks superposed on their wide, wind-broadened (optical) emission lines. Similar patterns of emission-line 
profile variations are detected in the Of stars ( Puppis and HD93129A, in the more evolved hydrogen-rich, 
luminous, Of-like WN stars HD93131 and HD93162, and in the more mass-depleted WC star in 7 2 Velorum. 
These observations strongly suggest that stochastic wind clumping is a universal phenomenon in the radiation- 
driven, hot winds from all massive stars, with similar clumping factors in all stages of mass depletion. 

Subject headings: line:profiles - stars: emission-line, Be - stars: mass-loss - stars: winds, outflows - stars: 
Wolf-Rayet 
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1. INTRODUCTION 

Evidence has been accumulating in the last two decades that 
hot stellar winds are far from being the smooth flows of es- 
caping matter that were often conveniently assumed. Cur- 
rent observations are of sufficient quality to show that the 
strong, hot stellar winds of bright, population I Wolf-Rayet 
(WR) stars are systematically pervaded by inhomogeneities 
on different scales. Variability on time scales of minutes to 
hours is so far best descr ibed by a fully clumped wind with 
no sm ooth component (Lepine & Moffatl [T999l : ILepine et al.\ 
1999). This clumping may very well be related to small-scale 
radiative instabilities lea ding to strong shocks in the wind 
(iGayley & Owockilll995h . 

Is this phenomenon unique to the extremely dense winds 
of WR stars, or do other hot stars also have similar globally 
clumped winds which have been difficult to detect due to their 
weaker emission lines? Subsequent to our WR studies, we 
found that the Of star ( Puppis shows ident ical structure in 
its wi nd, as found (statistically) in WR stars (Eversberg et al. 
1998): C Pup's optical He II A 4686 emission line varies in 
exactly the same manner as the isolated Hell 5411 line in 
WN stars or CM 5696 in WC stars. The variable emission 
component observed at very high S/N in £ Pup is also the 
same fraction of the total emission-line flux as in WR spectra, 
whose absolute amplitude is significantly less in O stars. This 
may explain why it had remained virtually unnoticed before 
in O stars. If this one O star is typical, then we expect to find 
the same stochastic emission-line variability in other O-star 
winds. 

Other types of hot stars also exhibit clumpy structure in 
their w inds. Examples includ e the [ WCl nuclei of plane tary 
nuclei (Grosdi dier et alj2003 ), nov ae (Lepi ne et af.Fl 999). su- 
pernovae (Mat hewson et al.\ 20001) and possibly others that 
have not been examined properly yet. 

Turbulent clumping in hot-star winds has a number of very 
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important consequences, e.g. among the most important are: 
(1) Most current estimates of mass-loss rates, which are sen- 
sitive to the square of the density, have to be revised down- 
wards by about a factor three with impo rtant implications for 
stellar evolution ilMoffat & Robert! [19911) . (2) Since clumps 
follow the mean flow, on average, they are turning out to 
be extremely useful in constraining the empirical wind ve- 
locity law, v(r), that is quite fundamental, but poorly known 
at present in hot-star winds; in particular, the convenient, 
popular inte rpolation /3-law appears to be only a rough ap- 
proximation (ILepine & Moffatll 19991) . (3) Wind clumps may 
ultimately provide the necessary means of compression and 
shielding required to form dust in some WR and nova winds 
(iWilhams et all [19871; ICherchneff et al.\ [19951 iBodel fl99l 
Marchenko & Moffat 2007). This might be relevant for dust 
formation in any type of star. The presumed associated de- 
crease in the mass-loss rates could have a major effect on their 
evolution, slowing down the shedding of their envelope and 
delaying their eventual transfer to a WR stage. 

One way out, however, is via the still poorly-understood, in- 
termediate, luminous blue variable (LBV) stage, which may 
provide the necessary strong episodic mass-loss to ultimately 
allow O stars to pass to WR stars dSmithl [2007b . In the 
meantime, massive-star models have been vastly improved. 
Including rotation in the models, which is generally quite 
rapid for O stars, will tend to feed more H-rich fuel into the 
core region and increase the luminosity for a given mass, 
thereby allowing more O stars to become WR, and relax- 
ing the requirement for unexplainably hi gh mass-loss rates 
in the HRD isochrone fits to observations (Maeder & Mevnet 
2000). Therefore, it appears particularly urgent to examine 
the clumping question in O-star winds of all types, so that evo- 
lutionary models for massive stars on/near the main sequence 
can be checked/constrained. This will increase our confidence 
in the age determinations of clusters and starbursts, which are 
turning out to be crucial tracers of star formation in the early 
Universe. It is also relevant in the question of the progenitors 
of the slow gam ma-ray bursts, believe d to arise in rapidly ro- 
tating WR stars (I van den Heuvelll2~007l) . The rotation rate will 
be higher if the progenitor loses less mass on the way to be- 
coming a WR star. In this context, it is interesting to note that 
clumping is still as strong in SMC WR stars at Zq /5 as at Zq 
in the Galaxy dMarchenko et al.l[2007l) . 

Very recent studies using the indirect approach have re- 
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TABLE 1 
Stars observed 



HD number 


Other name 


RA(2000) 


Dec(2000) 


V magnitude 


Spectral type 


Comment 


HD 66811 


£ Puppis 


08 03 35.04 


-40 00 11.3 


2.21 


04I(n)f 


single runaway 


HD93129A 




10 43 57.46 


-59 32 51.3 


6.97 


03If* 


visual binary 1 


HD 93131 




10 43 52.25 


-60 07 04.0 


6.48 


WN6ha 


single 


HD 93162 


7 2 Velorum 


10 44 10.33 


-59 43 11.4 


8.11 


WN6ha 


spectroscopic binary 2 


HD 68273 


08 09 31.95 


-47 20 11.7 


1.81 


WC8+07.5III-V 


spectroscopic binary 3 



Companion 0.05"from the primary i Nelan el al. 2004). 2 Orbital period 207.7 days i Gamen el al. 2006). 3 Orbital period 78.5 days i De Marco el al. 2000). 



vealed dec reased mass-loss rat es in O-star win ds by factors of 
3 or m ore dBouret et a/.ll2005l) to 10 or more (Fu llerton et al.\ 
2006). These are based on weak resonance and other lines in 
the far UV with the FUSE satellite. In the UV, the emission 
lines are stronger, but often saturated and in any case, it is dif- 
ficult to obtain sufficiently high S/N. A recent complementary 
approach in the UV has been successful, however, in reveal- 
ing variable features in the strong P Cygni absorption edges 
in the FUSE UV s pectru m of the single WC8 star WR135 by 
iMarchenko et al.\ (120061) . linking the presence of overdense 
clumps with observed shocks. 

In this study, we look for direct evidence for wind clump- 
ing in optical emission-lines generated in O-star and related 
winds. Clumping is detected in the acceleration zone near 
the base of the wind as transient, moving subpeaks on the 
emission-line profiles. Among O stars, Of stars are the easi- 
est to examine, since they have emission lines in the optical, 
where it is possible to obtain high S/N spectra of these nor- 
mally weak lines. 

2. OBSERVATIONS 

We observed two Of (one a repeat of previous observations 
of £ Pup) and two late-type, H-rich, weak-line, Galactic WN 
(super Of-like) stars intensively during three nights, 2000 Jan- 
uary 25/26-27/28, at the 4m telescope at CTIO. Apart from ( 
Pup, the other three stars are all located in the Carina Neb- 
ula. Parallel observations of the bright WR star 7 2 Vel (WC8 
+ Q7.5I II-V) were also ob served as a check star for known 
clumps (Lepine et al. 1999). Basic data on our five targets are 
found in Table 1 . 

We used the bench-mounted echelle spectrograph giving an 
inverse dispersion of 0.06 A per pixel and (2.7 -pixel) resolv- 
ing power R -30,000 over the AA 4200 - 5700 A range. The 
slit was 1" wide. Typical elementary exposures of 10-1000 
sec, depending on the stellar brightness, were combined from 
the CCD detector into one spectrum every —1-2 hours for 
each star, yielding net S/N « 500 per pixel (no binning) af- 
ter wavelet filtering (see below). Wavelengths were calibrated 
using Th-Ar spectra. An internal quartz-lamp was used for 
flat-fielding. 

The bright echelle B3V flux standard r\ Hya (Hainuv et al.\ 
fl992h was observed in an attempt to calibrate and splice to- 
gether the echelle orders; however, its standard flux values 
turned out to be of too low resolution to be useful, so its use 
was abandoned. This forced us to treat each echelle order sep- 
arately. In any case, this was not a serious impediment, as we 
chose to focus on the strongest, most isolated optical emission 
lines of Hell 4686 and Hell/H/J 4860 in each star. After some 
experimentation, we obtained final extracted spectra of good 
quality by dividing the spectra in order 9 by those in order 8 
for Hell 4686 and order 7 by order 6 for HeII/H/3 4860. In all 
the stars, orders 8 and 6 were essentially devoid of emission 



lines and thus could be used as continuum reference spectra. 
This division by neighboring orders allowed us to eliminate 
the strong order sensitivity-function, and rectify the spectra 
without polluting the target emission lines significantly. 

Nightly mean spectra for each of our 5 targets, are displayed 
in Figure 1 for the echelle orders covering the Hell 4686 and 
HeII/H/3 4860 atomic lines. All stars show those atomic lines 
strongly in emission, with large equivalent widths and very 
broad profiles which are typical of Wolf-Rayet and Of stars 
with strong, radiation-driven winds. Profiles from the Wolf- 
Rayet stars are dominated by the emission component, while 
profiles from the two Of stars show a slightly blue-shifted dip 
which is due to P Cygni type absorption from the stellar wind. 

3. ANALYSIS AND RESULTS 

We first examine whether there are significant global line 
variations from night to night in the P C ygni component, as 
seen before in the nightly s pectra of £ Pup (Moffat & Michaudl 
U98UlEversberg et aZ.ll998t) . The nightly means from each of 
the three nights are shown in Fig.l in different colors. Sig- 
nificant night-to-night variations are observed in the WC star 
7 2 Vel and in the two Of stars CPup and HD 93129A. Varia- 
tions are also apparent in the Hell 4686 line of the two WNha 
stars HD 93131 and HD 93162. The Hell 4686 line in C Pup 
is found to display the most pronounced variations, which 
suggest that the intens e line-profile variations discovered by 
Eversberg et al.\ (119981) might be an exceptional case. 

Variations on timescales of hours turn out to be much more 
interesting in all the stars. For each spectrum collected during 
the observing run, we subtracted the mean profile for the run. 
Despite all the calibration steps described in §2, many instru- 
mental artifacts remain in the spectra at the <1% level. While 
these artifacts have little effect on the mean spectra, they do 
show up prominently in the residuals. Most artifacts are ei- 
ther of low or high frequency in wavelength/pixel space. We 
thus filtered all of these out below 0.5 A and abov e —8 A us- 
ing wa velet transforms following the method of Lepi ne et al.\ 
( 1999). This is fortunate, given that most of the intrinsic vari- 
ations in Wolf-Rayet and Of star lines occur as narrow emis- 
sion line features (subpeaks) on a 1-2 A scale. 

In Figs. 2, 3, and 4, variability in the emission lines of Hell 
4686 (left) and HeII/H/3 4860 (right) are shown for each of 
our target stars. The mean profile for each line is displayed 
in the bottom panel. The middle panel plots the dispersion at 
each pixel for the entire observing run. The top panel is a two 
dimensional plot of the residuals (i.e. the instantaneous devi- 
ation from the mean profile) as a function of time and wave- 
length. Time runs upward on the vertical axis. The horizontal 
axis is labeled in units of radial velocity, calculated based on 
the redshift/blueshift of the emission line under consideration. 

Significant variations are observed in all the lines exam- 
ined. With variability occurring on timescales of hours, it is 
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FIG. 1. — Nightly mean spectra of our five program stars, listed in Table 1. For each star, two orders from the echelle spectrogram are shown, with the 
wavelength regions around the broad, intense emission lines of Hell 4686 (left) and HeII/H/3 4860 (right). The mean spectra for the nights of 2000 Jan 25, 26 
and 27 are shown in red, green and blue, respectively. 

clear that at least some of this variation will average out in the 
nightly mean. The night-to-night variations noted in Figure 
1 thus underestimate the variability in the line profiles. Con- 
versely, the mean from all three nights can be confidently used 
as a base profile, and the residuals from this means (as shown 
in Figs. 2, 3, and 4) will provide a good representation of the 
variable elements, be they either excess emission subpeaks or 
excess absorption troughs. 

In all the target stars, both the Hell 4686 and HeII/H/3 4860 
lines show variability which is extremely similar to that pre- 
viously observed in WR spectra. The same general pattern of 
propagation outwards from the line center is apparent in all 
the stars for at least one of the nights. These indicate out- 
wards acceleration in the winds of all these stars. Accordin g 
to the phenomenological models of lLepine & Moffal dl999). 
emission subpeaks in the line profile do not necessarily trace 
the acceleration of single, overdensities (clumps) in the wind. 
Rather, the emission subpeaks are likely the combined signa- 
ture of numerous clumps, so the outwards motion of subpeak 
elements across the line may not be a measure of the true wind 
acceleration. In any case, in the one star where the accelera- 
tion is most apparent (HD 93131) subpeak profiles on the blue 
and red edges are observed to move by about 400 km/s over 
the course of one night (8 hours). This suggests an outward 
wind-acceleration rate on the order of 14 m s" 2 , which is con- 
sistent ^wjfh_acc^krafenrates measured in other Wolf-Rayet 
stars dLepine & Moffatll 19991) . 

In Figure 5, we condense the variable information into plots 
of standard deviation per pixel across each line profile for all 
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FIG . 2 . — Time series of the residuals after subtraction of the mean profile, 
for the Wolf-Rayet star 7- Velorum. 
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Fig. 4.— 



Time series of the residuals after subtraction of the mean profile, for the Of stars f Puppis (left) and HD93129A (right). 




FIG. 5. — Comparison between the normalized mean profiles of the emission lines and the observed dispersion in the line profile, as a function of wavelength. 
Mean (thick) lines and dispersion (thin) lines are superposed on each plot, with the amplitude of the dispersion multiplied by a factor of 35 (hence a dispersion 
level of 0.35 on the plot actually corresponds to a dispersion amplitude of 0.01 in normalized flux units). 
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spectra collected over the three nights. Standard deviation 
profiles are all multiplied by a factor of 35, and compared 
with the mean profiles for each star. On this scale, the de- 
viation profiles and mean emission-line profiles have similar 
amplitudes, which suggest that all lines show variability at the 
r^3% level. This le vel of variability is on a par with previous 
data for 9 WR stars dLepine & Moffat! 19991). the Of star C Pup 
dEversbergef"^Hll998lh and ^ Vel dLepine et a/.lll999l) . The 
only exception is the He II A4860 profile in HD93162, whose 
variability is below the 1% level. 

4. DISCUSSION AND CONCLUSIONS 

Our spectroscopic time-series show systematic variability 
on a timescale of hours in bright emission-line profiles from 
hot, massive stars with moderately strong winds. The emis- 
sion lines examined here include lines which are significantly 
fainter, relative to the continuum, than the bright lines from 
Wolf-Rayet stars, where the phenomenon was first examined 
dLepine & Moffat! fl999l) . We find that the amplitude of the 
variations are at the l%-3% level of the line emission. This 
suggests that variability is universal in all emission lines, but 
that it is only apparent in the strongest of the lines. 

While the variability profiles (cr(A)), after rescaling, gen- 
erally follow the underlying line profiles (/(A)), they do not 
match exactly on close examination. There are two possi- 
ble causes for this: (1) the variability is strictly proportional 
to line emission, but the profile is distorted by an underlying 
absorption component, and/or (2) there exists an additional 
source of variability in the line profile, e.g. because of vari- 
ations in an under lying absorption c omponent, as in P Cygni 
absorption edges de.g. Robert 111992b . The effects from dis- 
torted line profiles are most apparent in the lines of the Of 
stars, which are flanked by a blueward absorption component, 
and are strongly asymmetric relative to the rest wavelength. 
Despite this, however, we clearly observe (Fig. 4) that the vari- 
ability profile is roughly symmetric about the rest wavelength, 
which confirms our assertion that the variability is associated 
mainly with the emission component. 

Overall, it appears that the clumps in the winds that cause 
the variable subpeaks across the emission lines, do indeed 
trace the wind. This would suggest that, within their turbu- 
lent speeds, the clumps obey essentially the same accelera- 
tion and velocity laws as the average wind itself. Indeed, the 
clump s are the wind! In fact, u sing extremely simple assump- 
tions, iLepine & Moffatl d 19991) were able to simulate strong 
WR emission lines and their observed variability solely us- 
ing a superposition of a very large number of discrete wind- 
emitting elements, without recourse to any kind of continuous 
background wind. This may not be proof of concept, but at 
least it is consistent with other astroph ysically turbulent me- 
dia such as the ISM ( e.g. StutzkHll993l) . We would thus claim 
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that hot-star winds are the (supersonic compressible) turbu- 
lent result of some kind of energy-input driving. Whether this 
occurs predominantly as a process of cascadin g dissipation 
from large to small scales (le.g. Henriksen|[T994T) or of shock- 
merging a ggregation fro m small radiati ve instabilities on sub- 
Sobolev dOwockill 19941) to large scales (Fel dmeier II 1 9951) . re- 
mains to be demonstrated observationally. 

The ratio er(A)//(A) is observed to vary between different 
lines, though only by a factor of a few at most for the same 
line in different stellar winds. The more pronounced case here 
is the star HD 193 162 in which Hell 4686 is more variable 
than HeII/H/3 4860 by about a factor of 3. This could be re- 
lated to enhanced Hell 4686 emission from the highly turbu- 
lent zone of wind-wind collision between the components in 
this massive binary system. The variability levels are how- 
ever remarkably similar between the various stars observed 
here. Since the intensity variations across optical (recombi- 
nation) spectral lines vary with the square of the density, the 
clumping filling factor, which varies as the square root of this 
quantity, must actually show much less dispersion compared 
to the observed line-intensity variations, from one wind to an- 
other. From this, we deduce that most hot-star winds (O, WR, 
LBV) require nearly the same reduction factor in mass-loss 
rate of typically ^3 (2-5 in the extreme) when the mass-loss 
rates are based on density-squared mechanisms. This result is 
in accord with the consensus on the ubiquity of the clumping 
phenomenon that came out of the rec ent workshop devote d 
entirely to clumping in hot-star winds (Hamann et al. 2008). 

Our conjecture is that spectroscopic observations with even 
higher signal-to-noise ratio should reveal similar variability 
patterns in all emission-line profiles of stars with hot winds. 
Wind clumping is the probable source of the variability. The 
timescale is set by the time it takes for individual clumps 
to propagate through the layer in the wind from which line 
emission occurs. Different lines will show different patterns, 
depending on the height and depth of the layer where emis- 
sion occurs. Emission lines formed in overlapping layers will 
have similar line-profile variability patterns. From a typical 
timescale of ^6 hours and a wind velocity speed of ~ 1000km 
s" 1 , the typical depth of the emission layers is ~2 10 7 km, or 
about 30 Rq, consistent with the idea that the lines are formed 
relatively close to the hydrostatic surface. All of these things 
should be tested with better future data. 
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